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^ . Abstract 

Os , 

^^ I We develop a general theory dealing with stochastic models for 

dynamical systems that are governed by various nonlinear, ordinary or 
partial differential, equations. In particular, we address the problem 
how flows in the random medium (related to driving velocity fields 
I^ ' which are generically bound to obey suitable local conservation laws) 

,^ . can be reconciled with the notion of dispersion due to a Markovian 

diffusion process. 



H Contexts 



Probabilistic concepts are ubiquitous in diverse areas of nonlinear science. 
Deterministic dynamical systems may give rise to random transport that 
is an intrinsic feature of their complex behaviour, augmented by a possi- 
ble choice of random initial or boundary data and/or suitable scaling limits. 
Nonequlibrium statistical physics directly employs random processes, Gaus- 
sian and non-Gaussian, which ultimately implement a nonlinear transport. 
Apart from a concrete identification of genuine sources of randomness in the 
dynamics of classical systems or in quantum theory, quite often we need 
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quantitative methods to deal with random-looking phenomena. Basically 
that refers to situations when origins of randomness are either uncontrollable 
(allowing merely for probabilistic predictions about the future behaviour of 
the system) or not definitely settled. 

The so-called Schrodinger boundary data and stochastic interpolation 
problem sets a conceptual and formal basis for a surprisingly rich group of 
topics. Here, IQ, stochastic analysis methods are used to deduce the most 
likely (generally approximate) underlying dynamics from the given (possibly 
phenomenological) input-output statistics data, pertaining to a certain dy- 
namical process that is bound to take place in a finite time interval. The 
pertinent motion scenarios range from processes arising in nonequilibrium 
phenomena, 0, ^, through classical dynamics of complex systems (deter- 
ministic chaos in terms of densities) to searches for a stochastic counterpart 
of quantum phenomena. They involve random processes that go beyond the 
standard Gaussian basis and enable a consistent usage of the jump- type pro- 
cesses (Levy and their perturbed versions) associated with the anomalous 
transport, 0]. In the diffusion processes context, we have identified the third 
Newton law for mean velocity fields as being capable to generate anomalous 
(enhanced) or non- dispersive diffusion- type processes through "perturbations 
of noise", 0. 

Since the stochastic interpretation of various differential equations is our 
major target, let us mention typical examples that are amenable to our 
methodology. Those are: Boltzmann, Navier-Stokes, Euler, Burgers (more 
or less standard transport of matter), Hamilton- Jacobi, Hamilton- Jacobi- 
Bellmann, Kardar-Parisi- Zhang (an issue of viscosity solutions and the in- 
terface profile growth), Fokker-Planck, Kramers (standard random propaga- 
tion related to the Brownian motion), both linear and nonlinear Schrodinger 
equations (probabilistic interpretation of solutions, in Euclidean and non- 
Euclidean versions of the problem, also with reference to an escape over a 
barrier and decay of a metastable state). 

In most of the above cases a natural linearisation of a nonlinear prob- 
lem is provided by generalized diffusion (heat) equations in their forward 
and/or backward versions, [|l|. On the contrary, a suitable coupled pair of 
time-adjoint nonlinear diffusion equations admits a linearisation in terms of 
the familiar Schrodinger equation. That involves Markovian diffusion pro- 
cesses with a feedback (enhancement mechanism named " the Brownian recoil 
principle"), ||3]. 



Concepts: Schrodinger' s interpolation prob- 
lem 

There are many procedures to deduce an intrinsic dynamics of a physical 
system from observable data. None of them is free of arbitrary assumptions 
needed to reduce the level of ambiguity and so none can yield a clean choice 
of the modelling strategy. As a standard example one may invoke the time 
series analysis that is a respected tool in the study of complex signals and 
is routinely utilised for a discrimination between deterministic and random 
inputs. 

Our objective is to reconstruct a random dynamics that is consistent 
with the given input-output statistics data. We shall outline an algorithm 
allowing to reproduce an admissible microscopic motion scenario under an 
additional assumption that the sought for dynamics actually is a Markovian 
diffusion process. This reconstruction method is based on solving the so- 
called Schrodinger boundary-data and interpolation problem. 

Given two strictly positive (usually on an open space-interval) boundary 
probability densities po{x), pt{x) for a process with the time of duration 
T > 0. One can single out a unique Markovian diffusion process which is 
specified by solving the Schrodinger boundary data problem: 

mT{A,B) = I d^x d^ymT{x,y) (1) 

Ja Jb 

^yn.ri3,V) = P.iS),j^.,nr(l,y)=PrW (2) 

where the joint probability distribution has a density 

rriT^x, y) = uoix) k{x, 0, y, T) vriy) (3) 

and the two unknown functions Uq{x), Vt{jJ) come out as (unique) solutions, 
of the same sign, of the integral identities. To this end, we need to have at our 
disposal a continuous bounded strictly positive (ways to relax this assumption 
are known) function k{x,s,y,t),0 < s < t < T, which for our purposes (an 
obvious way to secure the Markov property) is chosen to be represented by 
familiar Feynman-Kac integral kernels of contractive dynamical semigroup 
operators: 



k{y,s,x,t) = exp[- c((^(r), r)dr]o?/i}|''{(a^) (4) 



In the above, diJ,%'ll{uj) is the conditional Wiener measure over sample 
paths of the standard Brownian motion. (Another choice of the measure 
allows to extend the framework to jump-type processes.) 

The pertinent (interpolating) Markovian process can be ultimately deter- 
mined by means of positive solutions (it is desirable to have them bounded) 
of the adjoint pair of generalised heat equations: 

dtu{x, t) = u/\u{x, t) — c(x, t)u{x, t) (5) 

dtv{x,t) = —u/\v{x,t) + c{x,t)v{x^t) . (6) 

Here, a function c{x, t) is restricted only by the positivity and continuity 
demand for the kernel. 

Solutions, upon suitable normalisation give rise to the Markovian diffu- 
sion process with the factorised probability density p{x,t) = u{x,t)v{x,t) 
which, while evolving in time, interpolates between the boundary density 
data p{x,0) and p{x,T). The interpolation admits an Ito realisation with 
the respective forward and backward drifts defined as follows: 

r /^ X ^ Vu(x,t) _, 

U{X,t) 

in the prescribed time interval [0,T]. 

For the forward interpolation, the familiar Fokker-Planck (second Kolmogorov) 

equation holds true: 

dtpix, t) = uApix, t) - W\b{x, t)p{x, t)] (9) 

with p{x, 0) given, while for the backward interpolation (starting from p{x, T)) 
we have: 

dtp{x, t) = -uAp{x, t) - V[K{x, t)p{x, t)] . (10) 

The drifts are gradient fields, curlh = 0. As a consequence, those that 
are allowed by any prescribed choice of the function c(a:, t) must fulfill the 
compatibility condition 

c(f,t)=9t$ + ^(^ + V6) (11) 



which estabhshes the Girsanov-type connection of the forward drift b{x, t) = 
2z/V$(x, t) with the Feynman-Kac potential c{x, t). In the considered Schrodinger's 
interpolation framework, the forward and backward drift fields are connected 
by the identity b^, = b — 2z/V/np. 

For Markovian diffusion processes the notion of the backward transition 
probability density p*{y, s, x, t) can be consistently introduced on each finite 
time interval, say < s < t < T: 

p(f , t)p^{y, s, f , t) = p{y, s, X, t)p{y, s) (12) 

so that J p{y, s)p{y, s,x,t)d^y = p{x,t) and p{y, s) = J p*{y, s,x,t)p{x,t)d^x. 

The transport (density evolution) equations refer to processes running in 
opposite directions in a fixed, common for both time-duration period. The 
forward one executes an interpolation from the Borel set A to B, while the 
backward one executes an interpolation from B to A. 

Let us mention at this point that various partial differential equations 
associated with Markovian diffusion processes are known not to be invariant 
under time reversal. That implies dissipation and links them with irreversible 
physical phenomena. However, the correspoding processes are known to ad- 
mit a statistical inversion and asking for a statistical past of the process 
makes sense. 

In particular, the knowledge of the Feynman-Kac kernel implies that the 
transition probability density of the forward process reads: 

p{y,s,x,t) = k{y,s,x,t) .J . ■ (13) 

v{y,s) 

while the corresponding transition probability density of the backward pro- 
cess has the form: 

lliv S) 

P* [y, s, X, t) = k{y, s, X, t) J . (14) 

U[X,t) 

Obviously in the time interval < s < t <T there holds: 

u{x, t) = J uo{y)k{y, s, x, t)d^y (15) 

v{y,s)= k{y,s,x,T)vT{x)d^x . (16) 



Consequently, we have fully determined the underlying (Markovian) ran- 
dom motions, forward and backward, respectively. All that accounts for 
perturbations of (and conditioning upon) the Wiener noise. 



Particularities 

If we consider a fluid in thermal equilibrium as the noise carrier, a kinetic 
theory viewpoint amounts to visualizing the constituent molecules that col- 
lide not only with each other but also with the tagged (colloidal) particle, so 
enforcing and maintaining its observed erratic motion. The Smoluchowski 
approximation takes us away from those kinetic theory intuitions by project- 
ing the phase-space theory of random motions into its configuration space 
image which is a spatial Markovian diffusion process, whose formal infinites- 
imal encoding reads: 

— * 

dX(t) = -^dt + V2DdW(t) . (17) 

mp 

In the above m stands for the mass of a diffusing particle, /9 is a friction 
parameter, D is a diffusion constant and W{t) is a normalised Wiener process. 
The Smoluchowski forward drift can be traced back to a presumed selective 
action of the external force F = —W on the Brownian particle that has 
a negligible effect on the thermal bath but in view of frictional resistance 
imparts to a particle the mean velocity F /mj3 on the /3~"^ time scale. 

The noise carrier (fluid in the present considerations) statistically remains 
in the state of rest, with no intrinsic mean flows, hence unperturbed in the 
mean (all other cases we associate with the term "perturbations of noise"). 
At the first glance, a formal replacement of the Smoluchowski forward drift 
in Eq. (17) by any space-time dependent driving velocity field would suggest 
a legitimate procedure to implement a net transport that would combine 
dispersion due to a diffusion process with deterministic mean flows due to 
external agencies (consider Euler, Navier-Stokes or Burgers velocity flelds as 
reference examples). However the situation is not that simple. 

It is well known that a spatial diffusion (Smoluchowski) approximation 
of the phase-space process, allows to reduce the number of independent local 



conservation laws to two only. Therefore the Fokker-Planck equation can al- 
ways be supplemented by another (independent) partial differential equation 
to form a closed system. If we assign a probability density Pq{x) with which 
the initial data xq = X{0) are distributed, then the emergent Fick law would 
reveal a statistical tendency of particles to flow away from higher probabil- 
ity residence areas. This feature is encoded in the corresponding Fokker- 
Planck equation: dtp = —V ■ (vp) = —V ■ [{:£g — D^)p] where a diffusion 

current velocity is v{x,t) = b{x,t) — D '^S^'J while the forward drift reads 

b{x,t) = -^. Clearly, the local diffusion current (a local flow that might be 

experimentally observed for a cloud of suspended particles in a liquid) j = vp 
gives rise to a non-negligible matter transport on the ensemble average, even 
if no intrinsic mean flows are in existence in the random medium. 

We may surely consider a formal replacement of the previous ^ by a 

local velocity field v{x,t). However, irrespectively of whether we utilize ^ 
or v{x, t) in the formalism, the velocity field must obey the natural (local) 
momentum conservation law which directly originates from the rules of the 
Ito calculus for Markovian diffusion processes and from the first moment 
equation in the diffusion approximation (!) of the Kramers theory: 

dtV + {v-V)v = V{yi-Q) . (18) 

An effective potential function VL{x) can be expressed in terms of the 
Smoluchowski forward drift bix) = — ^ as follows: il = TT^^rm + -^V ■ F. 
That is to be compared with Eq. (11) which governs the more general space- 
time dependent situation. 

Moreover, in the would-be Euler equation (18), instead of the standard 
pressure term, there appears a contribution from a probability density p- 
dependent potential Q{x,t). It is given in terms of the so-called osmotic 
velocity field u{x,t) = DV lnp{x,t): Q{x,t) = \u^ + DV ■ u and is generic 
to a local momentum conservation law respected by isothermal Markovian 
diffusion processes. 

To analyze general perturbations of the medium and then the resulting 
intrinsic (mean) flows, a function 6(X(t),t), must replace the Smoluchowski 
drift. Under suitable restrictions, we can relate probability measures corre- 
sponding to different (in terms of forward drifts !) Fokker-Planck equations 



and processes by means of the Cameron-Martin-Girsanov theory of measure 
transformations, [|l[. For suitable forward drifts which are gradient fields 
that yields the most general form of an auxiliary potential (cf. Eq. (11)) 
n{x, t) = 2D[9t0 + i(^ + V ■ b)] . We denote b{x, t) = 2DV<p{x, t). 

Mathematical features of the formalism appear to depend crucially on 
the properties (like continuity, local and global boundedness, Rellich class) 
of the auxiliary potential Vt. Let us consider a bounded from below (lo- 
cal boundedness from above is useful as well), continuous function f2(x, t). 
Then, by means of the gradient field ansatz for the diffusion current velocity 
{v = VS -^ dtp = —V ■ [(V5')p]) we can transform the local momentum con- 
servation law of a Markovian diffusion process to a universal Hamilton- Jacobi 
form: 

n = dtS+^\WS\^ + Q (19) 

where Q{x,t) was defined before. By applying the gradient operation we 
recover the conservation law (18). 

In the above, the contribution due to Q is a direct consequence of an initial 
probability measure choice for the diffusion process, while Q does account for 
an appropriate forward drift of the process via Eq. (19). 

Thus, in the context of Markovian diffusion processes, we can consider a 
dosed system of partial differential equations which comprises the continu- 
ity equation dtp = —V{vp) and the Hamilton- Jacobi equation, plus suitable 
initial (and/or boundary) data. The underlying isothermal diffusion process 
is specified uniquely. Those two partial differential equations set ultimate 
limitations in the old-fashioned problem of "how much nonlinear" and "how 
much space-time dependent" the driving velocity field can be to yield a con- 
sistent stochastic diffusion process (or the Langevin-type dynamics) 
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